Kneelike structure in the spectrum of the heavy component of cosmic rays observed 

with KASCADE-Grande 
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We report the observation of a steepening in the cosmic ray energy spectrum of heavy primary 
particles at about 8 ■ 10^® eV. This structure is also seen in the all-particle energy spectrum, but is 
less significant. Whereas the 'knee' of the cosmic ray spectrum at 3-5 ■ 10^^ eV was assigned to light 
primary masses by the KASCADE experiment, the new structure found by the KASCADE-Grande 
experiment is caused by heavy primaries. The result is obtained by independent measurements of 
the charged particle and muon components of the secondary particles of extensive air showers in the 
primary energy range of 10^® to lO'^^ eV. The data are analyzed on a single-event basis taking into 
account also the correlation of the two observables. 

PACS numbers: 98.70.Sa, 95.85.Ry, 96.50.Sd 
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The determination of the primary energy and com- 
position in the energy range from 10^^ eV up to above 
10^° eV has been subject of earthbound experiments for 
more than five decades. It has been shown that the 
high-energy all-particle spectrum has a power-law like 
behavior (oc E'^ , 7 w —2.7), with features known as the 
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'knee' at 3-5 • IQi^eV and the 'ankle' at 4-10 • IQi^eV, 
respectively. Whereas at the knee the spectrum steep- 
ens (A7 = —0.3-0.4), the ankle is characterized by a 
flattening of the spectrum (A7 = -^ 0.3-0.4). The KAS- 
CADE experiment has shown that the knee is due to a 
distinct decrease in the Hxix of primaries with light mass 
(Z < 6) pj, l2[. Many astrophysical models comment- 
ing on the origin of the knee assume a dependence of 
such break-offs ('knees') on the charge of the primary 
nuclei [3|, Ij]. Assuming that the knee is related to a 
break in the spectrum of primary Hydrogen nuclei, a 
knee-like structure in the spectrum of the heavy compo- 
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FIG. 1: Layout of the KASCADE-Grande experiment: shown 
are the Grande array as well as the KASCADE array with its 
central detector (CD) and muon tracking detector (MTD). 
The shaded area marks the outer 12 clusters (16 detector 
stations each) of the KASCADE array consisting of shielded 
(muon array) and unshielded detectors. The inner 4 clusters 
consist of unshielded detectors, only. 



nent (Z > 13 up to Iron nuclei) is expected in the energy 
range from about 4 ■ 10^^ eV to about 1.2 ■ lO^^eV. So 
far, such a structure has not been observed experimen- 
tally. We present measurements of extensive air showers 
(EAS) in the primary energy range of 10^^ eV to 10^* eV 
performed with KASCADE-Grande (KArlsruhe Shower 
Core and Array DEtector with Grande extension) and 
investigate the mass composition of the cosmic rays. 

KASCADE-Grande, located at 49.1°N, 8.4°E, 
llOma.s.L, consists of the Grande array with 37 stations 
of 10 m^ scintillation detectors each, spread over an area 
of 700 X 700 m^, the original KASCADE array covering 
200 X 200 m^ with unshielded and shielded detectors, 
a muon tracking device, and a large calorimeter [5|, |6|. 
This multi-detector system allows us to investigate in 
detail the EAS generated by high-energy primary cosmic 
rays in the atmosphere. For the present analysis, the 
estimation of energy and mass of the primary particles 
is based on the combined measurement of the charged 
particle component by the detector array of Grande and 
the muon component by the KASCADE muon array 
(Fig. [T]). Basic shower observables like the core position, 
zenith angle, and total number of charged particles 
(shower size Nch) are derived from the measurements 
of the Grande stations. While the Grande detectors are 
sensitive to charged particles, the muonic component is 
measured independently by the shielded detectors of the 
KASCADE array. 192 scintillation detectors of 3.24 m^ 
sensitive area each are placed below an iron and lead 
absorber to select muons above 230 MeV kinetic energy. 
A core position resolution of 5 m, a direction resolution 
of 0.7°. and a resolution of the shower size of about 15% 
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FIG. 2: Two-dimensional distribution of the shower sizes: 
charged particle number and total muon number. All qual- 
ity cuts are applied. In addition, a roughly estimated energy 
scale is indicated in the upper panel. The lower panel shows 
a zoom to higher energies of the same observables, now cor- 
rected for attenuation. 



are achieved. The total number of muons (A^^) with a 
resolution of about 25% is calculated by combining the 
core position determined by the Grande array and the 
muon densities measured at the KASCADE array, where 
Nf^ undergoes a correction for a bias in reconstruction 
due to the asymmetric position of the detectors [5| . 

The present analysis is based on 1173 days of data 
taking. The cuts on the sensitive area (EAS core re- 
constructed within the array) and zenith angle (< 40°), 
chosen to assure best and constant reconstruction accu- 
racies, result in an exposure of 2 • 10^'^ m^-s-sr. Figure [2] 
displays the correlation of the two observables Nch and 
Nfj^. This distribution is the basis of the following anal- 
ysis, since it contains all the experimental information 
required for reconstructing energy and mass of the cos- 
mic rays: the higher the energy of the primary cosmic 
ray the larger the total particle number. The fraction 
of muons of all charged particles at observation level is 
characteristic for the primary mass: showers induced by 
heavy primaries start earlier in the atmosphere and the 
higher nucleon number leads to a relatively larger muon 
content at observation level. KASCADE-Grande mea- 
sures the particle number at an atmospheric depth well 
beyond the shower maximum, where the electromagnetic 
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FIG. 3: Evolution of the k parameter as a function of the re- 
constructed energy for experimental data compared with sim- 
ulations of primary masses for the angular range 0-24° . The 
error bars assign statistical as well as reconstruction uncer- 
tainties of k. The line displays the chosen energy dependent 
fc-values for separating the mass groups, where the dashed 
lines assign the uncertainty of the selection. 



component already becomes reduced. Thus, electron-rich 
EAS are generated preferentially by light primary nuclei 
and electron-poor EAS by heavy nuclei, respectively. 

However, a straightforward analysis is hampered by 
the shower-to-shower fluctuations, i.e. by the dispersion 
of the muon and electromagnetic particle numbers for a 
fixed primary mass and energy. In addition, cosmic rays 
impinging on the atmosphere under different zenith an- 
gles show a varying, complicated behavior due to the non- 
uniform mass and density distribution of the air. There- 
fore, the absolute energy and mass scale have to be in- 
ferred from comparisons of the measurements with Monte 
Carlo simulations. This creates additional uncertainties, 
since the physics of the relevant particle interactions is 
not completely tested by man-made accelerator experi- 
ments. The uncertainties imposed by the hadronic inter- 
action models are more relevant for composition analy- 
ses than for energy measurements. Hence, our strategy 
is to separate the measured EAS in electron-poor and 
electron-rich events as representatives of the heavy and 
light primary mass groups, similar to the analysis pre- 
sented in Ref. [1|. The shape and structures of the re- 
sulting energy spectra of these individual mass groups 
are much less affected by the differences of the various 
hadronic interaction models than the relative abundance. 

As a consequence of the considerations above, the en- 
ergy and mass assignment of individual events is achieved 
by using both observables Nch and TV^, as well as their 
correlation. The following equation is motivated by dis- 
cussions of hadronic air showers in reference [7|, with 
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FIG. 4: Reconstructed energy spectrum of the electron-poor 
and electron-rich components together with the all-particle 
spectrum for the angular range 0-40°. The error bars show 
the statistical uncertainties; the bands assign systematic un- 
certainties due to the selection of the subsamples. Fits on the 
spectra and resulting slopes are also indicated. 



the basic idea that the total number of secondary par- 
ticles at observation level is related to the primary en- 
ergy while the energy sharing between the electromagne- 
tic and the hadronic (i.e. muonic) shower components is 
related to the primary mass. Therefore, the primary en- 
ergy log]^g(_B) is assumed to be proportional to the shower 
size \ogiQ{Nch) with a correction factor that accounts for 
the mass dependence by making use of the measured ra- 
tios of shower sizes \ogiQ{Nch/N^): 
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with logio(7Vcft/A^^)ff,Fe = CH,Fe ' logio{Nch) + dn^Fe- 

The parameter k takes into account both the average dif- 
ferences in the Nch/N^ ratio among different primaries 
with same Nch as well as the shower to shower fluctu- 
ations for events of the same primary mass. The exact 
form of the equation is optimized for the experimental sit- 
uation of KASCADE-Grande and the free parameters [8| 
are determined by Monte Carlo simulations [16[. They 
are defined independently for 5 different zenith angle in- 
tervals of equal exposure (the upper limits of 9 are 16.7°, 
24.0°, 29.9°, 35.1°, and 40.0°) to take into account the 
shower attenuation. Data are combined only at the very 
last stage to reconstruct the final energy spectrum. The 
-/Vch-A^/i-correlation of individual events is incorporated 
in calculating k, which serves now as mass sensitive ob- 
servable. Fig. [3] shows the evolution of fc as a function 



of the reconstructed energy for the first two zenith angle 
bins, where a similar behavior is observed for all angular 
ranges. The error bars include statistical as well as recon- 
struction uncertainties of the fc-parameter. The width of 
the k distributions decreases slightly for increasing en- 
ergy and amounts, at 100 PeV, to about ±0.2, ±0.15, 
±0.4 for H, Fe and data, respectively. 

The fc-parameter is used to separate the events into 
different samples. The line in Fig. |3] separates the 
electron-poor (heavy) group, and is defined by fitting the 
kepiE) = {ks^iE) + kc{E))/2 distribution. The dashed 
lines represent the uncertainties in defining this energy 
dependent selection cut. The resulting spectra are shown 
in Fig. m where the band indicates changes of the spec- 
tra when the cut is varied within the dashed lines shown 
in Fig. |3l The energy resolution for an individual event 
is better than 25% over the entire energy range and the 
all-particle spectrum is reconstructed within a total sys- 
tematic uncertainty in flux of 10-15% [8| [17| . 

The reconstructed spectrum of the electron-poor 
events shows a distinct knee-like feature at about 8 • 
10"'^^ eV. Applying a fit of two power laws to the spec- 
trum interconnected by a smooth knee [l3| results in 
a statistical significance of 3.5(7 that the entire spec- 
trum cannot be fitted with a single power law. The 
change of the spectral slope is A7 = —0.48 from 7 = 
-2.76 ± 0.02 to 7 = -3.24 ± 0.05 with the break po- 
sition at \og^Q{E/eV) = 16.92 ± 0.04. Applying the 
same function to the all-particle spectrum results in a 
statistical significance of only 2.1(t that a fit of two 
power laws is needed to describe the spectrum. Here 
the change of the spectral slope is from 7 = —2.95 ± 0.05 
to 7 = —3.24 ± 0.08, but with the break position again 
at logio(S/ey) = 16.92 ± 0.10. Hence, the selection of 
heavy primaries enhances the knee-like feature that is al- 
ready present in the all-particle spectrum. The spectrum 
of the electron-rich events (light and medium mass pri- 
maries) is compatible with a single power law with slope 
index 7 = —3.18 ± 0.01. However, a recovery to a harder 
spectrum at energies well above 10"'^'' eV cannot be ex- 
cluded by the present data. This finding is of particular 
interest and needs more detailed investigations with im- 
proved statistics in future. 

The main result, i.e. the knee- like structure in the spec- 
trum of electron-poor events, is validated in the follow- 
ing by various cross-checks (Fig. ^. Variations of the 
slopes of the selection cut, as well as parallel shifts of 
the cut lines have shown that the spectral form, i.e. the 
knee-like structure of the electron-poor event sample, is 
retained. By shifting k to larger values the fraction of 
heavy primaries in the sample is enriched. Interestingly, 
we found that the slope index of the spectrum is not sig- 
nificantly changing beyond the break, but gets systemat- 
ically harder at lower energies . The position of the break 
remains constant, indicating that the heaviest primaries 
in the sample dominate the spectral form. An example of 
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FIG. 5: Energy spectra of electron-poor (heavy) event sam- 
ples obtained by different selection and reconstruction crite- 
ria. The original spectrum from Fig. |4] is compared with the 
spectrum from a more selective cut in the fc-parameter and 
with the spectrum obtained by using the Vc/c-parameter for 
selecting the electron-poor events (see text). 



a spectrum obtained by such a variation of the selection 
cut is shown in Fig. [5l 

A systematic uncertainty possibly affecting the inter- 
pretation of the data is related to the attenuation of the 
particle numbers in the atmosphere. So far, the atten- 
uation given by the EAS simulations is taken into ac- 
count. For validation, an independent analysis is per- 
formed where the correction for attenuation, i.e. for the 
zenith angular dependence, is based on the measured 
events, and not on simulations. The correction parame- 
ters arc obtained by applying the Constant Intensity Cut 
Method (CIC) [lj| to the two obscrvablcs independently. 
This procedure allows the data collected from different 
zenith angles to be combined in a model independent way. 
The shower size ratio Ycic — log^Q N'^/ log^g iV^^ is cal- 
culated, where N' and N'^^^ are the shower sizes corrected 
for attenuation effects in the atmosphere in such a way 
that they correspond to the shower sizes at a certain ref- 
erence zenith angle. In order to check, in addition to the 
attenuation correction, also reconstruction and selection 
uncertainties, we applied for this analysis more stringent 
cuts, which increase the energy threshold and decrease 
the statistics of the event sample compared to the stan- 
dard analysis. Now, Ycic is used to separate the events 
into electron-rich and electron-poor subsamples. In con- 
trast to the fc-parameter, the Yc/c-parameter is almost 
energy independent, where the energy of the individual 
events is again determined using eqn. [TJ For direct com- 
parison with the results obtained before, Ycic > 0.845 
is chosen for selecting the electron-poor event sample. 
The reconstructed spectrum (see Fig. [5]) obviously con- 
firms the earlier finding of the knee-like structure, which 
is due to a decrease in the flux of the heavy component. 

Another source of systematic uncertainty is related to 



the hadronic interaction model. In the frame of QGSJet- 
II, the measured distributions in k and Ycic a-re in agree- 
ment with a dominant electron-poor composition for the 
entire energy range. Whereas the Ycic and k values 
themselves behave differently for other hadronic inter- 
action models, the measured and simulated Yqic- and k- 
dependences on energy, and hence the shapes and struc- 
tures of the resulting spectra, are similar [18J. Details 
will be discussed in a forthcoming paper, but it is not 
expected that the basic result of the present analysis 
changes. 

Summarizing, by dividing KASCADE-Grande mea- 
sured air-shower events in electron-rich and electron-poor 
subsamples, there is first evidence that at about 8-10^^ cV 
the spectrum of the heavy component of primary cosmic 
rays shows a knee-like break. The spectral steepening 
occurs at an energy where the charge dependent knee 
of primary iron is expected, when the knee at about 3- 
5 • 10"'^^ eV is assumed to be caused by a decrease in the 
flux of primary protons. 
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